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ABSTRACT 

We present a method for the computation of the variance of cosmic microwave back- 
ground (CMB) temperature maps on azimuthally symmetric patches using a fast con- 
volution approach. As an example of the application of the method, we show results 
for the search for concentric rings with unusual variance in the 7- year WMAP data. 
We re-analyse claims concerning the unusual variance profile of rings centred at two 
locations on the sky that have recently drawn special attention in the context of the 



conformal cyclic cosmology scenario proposed by [Penrose (2009). We extend this anal- 
ysis to rings with larger radii and centred on other points of the sky. Using the fast 
convolution technique enables us to perform this search with higher resolution and a 
wider range of radii than in previous studies. We show that for one of the two special 
points rings with radii larger than 10 degrees have systematically lower variance in 
comparison to the concordance ACDM model predictions. However, we show that this 
deviation is caused by the multipoles up to order i = 7. Therefore, the deficit of power 
for concentric rings with larger radii is yet another manifestation of the well-known 
anomalous CMB distribution on large angular scales. Furthermore, low variance rings 
can be easily found centred on other points in the sky. In addition, we show also the 
results of a search for extremely high variance rings. As for the low variance rings, 
some anomalies seem to be related to the anomalous distribution of the low order 
multipoles of the WMAP CMB maps. As such our results are not consistent with the 
conformal cyclic cosmology scenario. 

Key words: methods: data analysis — cosmic background radiation — cosmology: 
observations 



1 INTRODUCTION 

A wide range of cosmological projects carried out in the last 
two decades - observation of supernova, galaxy surveys and 
measurements of cosmic microwave background (CMB) fluc- 
tuations - have resulted in high quality data enabling the es- 
tablishment of a standard cosmological model i.e. the ACDM 
model. This model is based on the Big Bang paradigm for- 
mulated in the last century by Alp her, Her man &; Ga mow 
(1948} and extended by the inflationary paradigm by Guth 



( 1981| . It explains the observed large scale structure of the 



Universe as the result of the growth of primordial inhomo- 
geneities in an expanding Universe, generated during an in- 
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flationary epoch, and driven by the mechanism of gravita- 
tional instability. The rate of expansion and the growth of 
inhomogeneities depends on the relative amounts of the dif- 
ferent forms of matter and energy that constitute the total 
content of the Universe. In the present epoch the Universe 
consists of 5% of baryonic matter, 22% of dark matter and 
73% of dark energy ( |Larson et al.||2011| . With only a mod- 
est number of free parameters, the ACDM model is able to 
successfully fit thousands of observational data points. 



However, this model can not be considered as a funda- 
mental model. It does not explain the nature of dark matter 
nor dark energy. Furthermore, it postulates the existence of 
an unknown scalar field responsible for the inflation. Thus, 
there is a need for stringent tests of the validity of this model 
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by comparison with alternative cosmological theories having 
different observational predictions. 

One example of such an approach is the search in CMB 
maps for features that are not predicted by an isotropic 
and homogeneous standard cosmological model. Remark- 
ably, this effort has resulted in several reports of both non- 
Gaussianity and evidence for the breakdown of statisti- 
cal isotropy in the Wilkinson Microwave Anisotropy Probe 
( WMAP) CMB data, as established by many qualitatively 
different methods. In particular, an intriguing planarity and 
alignment of the quadrupole and octopole has been demon- 
strated (|de Oliveira-Costa et al.||2004| |Copi et al.||2004| 



Schwarz et al.|2004| |Bielewicz et al.|2005| |Land &; Magueij o 



2005). A localised source of non-Gaussianity, in the form of 
a very cold spot on the sky of angular scale 10 degrees, has 
been detected by|Vielva et al.|(|2QQ4]> andjCruz et al.|(|2005|). 



In addition, an asymmetry in large-scale power as measured 
in the two hemispheres of a reference frame close to that 



delineated by the ecliptic plane was determi ned by|Eriksen 
et al.| (|2004| and [Hansen et al.| ( |2004| (also |Hansen et al. 
Hoftuft et al 



2009; 



2009) 



An alternative approach is to search the CMB maps for 
specific patterns that are signatures of alternative cosmolog- 
ical models. One example of such signatures are the pairs of 
circles with matching distribution of the CMB anisotropy 
signal predicted for universes with multiconnected topol- 
ogy ( |Cornish, Spergel, fe Starkman||1998l ). If the size of the 
fundamental domain were smaller than diameter of the ob- 
servable Universe then we should see at least one pair of 
such matched circles. Such searches, as most recently un- 



dertaken by |Bielewicz &; Ba nday (2011), have yielded nega- 
tive results (see also Key et al.||2007 ). However, as pointed 
out by |Bielewicz, Banday, &; Gorski (2012), with the up- 



coming Planck data we should be able to perform a similar 
search using CMB polarisation data as a cross-check of the 
temperature-based results (see also |Riazuelo et al.| 2006). 

A second example is provided by models of pre- 
inflationary massive particles ( |Fialkov, Itzhaki, fc Kovetz 
|2010| ). These predict the existence of concentric cold and hot 
rings in the CMB and a bulk flow of galaxies toward their 
center. The detection of such rings for the WMAP maps at 
3cr confidence level was reported recently by Ko vetz, Ben-| 



David, & Itzhaki (2010) 



Another example of a model that imprints circular pat- 
terns on the CMB sky is due to the theory of eternal infla- 
tion (Fee ney et al.|2011 ). Here, collisions of bubble universes 
leave signatures in the CMB maps with a characteristic pro- 
file. Hints of evidence of such signatures in the WMAP data 



were recently reported by |Feeney et al.| ( |2011| and |McEwen| 



et al. 



(2012). 



Finally, the scenario of conformal cyclic cosmology pro- 
posed by |Penrose] ( 2009 ) (see also |Penrose||2008[ [2010 ) pro- 
poses that the collisions of massive black holes in a previous 
cycle of the Universe results in a set of concentric rings with 
low variance on the CMB sky. In this context, rings centred 
at points with the Galactic coordinates (l,b) = (105°, 37°) 
and (l,b) — (252°,— 31°) have recently merited special at- 
tention jHajian||2011| |Moss, Scott, fe Zibin|[2011| |Wehus~ 



Eriksen|201~ 

In this paper, we perform a complete search for rings 
with anomalously variances in the WMAP CMB maps. For 
the computation of variance, we use an approach based on 



fast convolution on the sphere thus enabling the search for 
concentric rings centred on a high resolution grid and with 
radii in the range from 0° to 90°. This significantly extends 
previous studies that were limited, due to the use of pixel 
based methods, to rings centred on a lower resolution grid 
and radii in the range [0°,20°]. Furthermore, in addition 
to the low variance rings we search also for the high vari- 
ance rings. Though, they are not predicted by the conformal 
cyclic cosmology scenario, their presence would manifestly 
contradict the standard inflationary models. 

It should be noted that |McEwen et al. ( 2012| used a 
similar approach in the context of searching for signatures 
of eternal inflation. 



2 DATA 

The search for extreme variance rings was performed on the 
7-year WMAP data ( |Jarosik et al.|2011| ). The WMAP satel- 
lite observes the sky in five frequency bands denoted K, Ka, 
Q, V and W, centred on the frequencies of 22.8, 33.0, 40.7, 
60.8, 93.5 GHz with angular resolutions of approximate! 
52.8', 39.6', 30.6', 21' and 13.2', resp ectively. The map:' 
are pixelised in the healpdQ scheme ( Gorski et al. 2005 ) 
with a resolution parameter N s id e = 512, corresponding to 
3,145,728 pixels with a pixel size of ^ 7 arcmin. In par- 
ticular, we studied the V and W-band maps corrected for 



Galactic emission using a template fitting approach (Gold 
eTaLl[20lT| . 



3 COMPUTATION OF VARIANCE BY FAST 
CONVOLUTION ON THE SPHERE 

The variance of the CMB in concentric rings with radius r 
and width Ar centred on pixel i, Var*, r (AT), can be esti- 
mated via the expression: 



Var*, r (AT) = 



J2j ATjMjBjj / Y^j ATjMjB, 

J2k M k B lk ' 



+ 



(i) 



where Mi denotes mask applied to the map and is the 
ring profile such that 



Bh = { J 



for r < arccos(n; • hj) < r + Ar 
otherwise 



(2) 



The sum in the denominators is over the unmasked pixels in 
a given ring. The last term in equation |l]) corresponds to a 
noise variance term that has to be subtracted to yield an un- 
biased estimation of the CMB variance. The noise variance 
in a given pixel cr 2 is defined as a 2 = o-l hs /N° hs , where a Q bs 
is the rms of noise per observation and A/"° bs denotes number 
of observation of a given pixel. Both of these quantities are 
provided by the WMAP team for each map. 

The correction for the noise bias is not needed if we 
estimate the cross-variance of two differencing assemblies 
(DAs) for a given band. Such an estimator is not susceptible 
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to uncertainties in the modeling of the noise properties. In 
this case, the estimator of the variance takes the form: 



Vari, r (AT) 



+ 



(3) 



>: / AV,.,A7, J ,\/,/r, 

(£. ay',.^/, /;;;,•) (v a />>;•) 



{E k M k Bi k y 

where ATi and AT2 denote maps for two DAs. We will use 
this estimator to derive most of the results presented below. 
We checked that the differences between the two estimators 
are negligible for the 7-year WMAP data. In the case of 
the W-band, we used coadded maps for two pairs of DAs 
i.e. coadded (Wl+W2)/2 and (W3+W4)/2 maps. 

The computation of the estimator of CMB variance by 
direct summation over the pixels within the ring for all ring 
centers is very time consuming. However, one can observe 
that the variance estimator given by equation is simply 
the convolution of the square of the masked map, corrected 
for the mean, with the ring profile of a given radius and 
width. For the cross- variance, the product of the two maps is 
convolved with the ring profile. Thus, the computations can 
be significantly speeded up by use of a technique using fast 
convolution on the sphere and therefore performing all essen- 
tial computations in the spherical harmonics space. Because 
the ring profile is azimuthally symmetric the computations 
are the same as for convolution with an azimuthally sym- 
metric beam. In order to estimate the variance, we need to 
perform four such convolutions. We used the fast spherical 
harmonics transform implemented in the HEALPIX scheme 
for these computations. 

Such an approach allows us to perform a search for rings 
with an unusual variance centred on a higher resolution grid 
and for a wider range of ring radii than by using the pixel 
based methods employed in previous studies. It can be also 
easily extended to higher moments of the map such as skew- 
ness and kurtosis, as well as used for any azimuthally sym- 
metric profiles such as discs with an arbitrary radial weight- 
ing function. 

Examples of the spherical harmonic window functions 
for example ring geometries, bg = Bio / \/2£ + 1 , are shown 
in Fig. ^ Due to the azimuthal symmetry of the ring, the 
spherical harmonics coefficients, 

rcos(r) 

Bim oc 5 m o / Pi(cos6) d cos # , (4) 

J cos(r+Ar) 

where Pi(cosQ) are the Legendre polynomials, are non-zero 
only for modes m = 0. The factor \/2£ + 1 comes from the 
averaging of the squares of the coefficients over m. The win- 
dow functions were normalised to one at 60. Note that, due 
to pixelisation effects, the ring window function correspond- 
ing to a pixelised sphere is not exactly equal to the analyti- 
cal expression for the window function given by relation Q . 
Therefore, in our computations we used the former since the 
data is pixelised in the healpix scheme. 

As can be seen, the ring window function oscillates 
around zero with a wavelength ~ 2n/r dependent on the 
ring radius r expressed in radians i.e. 36 and 9 for rings 
of width 0.5° with outer radii of 10° or 40°, respectively. 
Furthermore, because of the v / 2?-FT factor in the denomi- 
nator, the envelope of the widow function falls as £~ 1 ^ 2 . Ad- 
ditional modulation of the envelope depends on the width of 



the ring Ar. To better understand this behaviour, the ring 
profile can be considered as a disc of given outer radius from 
which a smaller disc of the given inner radius is subtracted. 
The interference of the window functions with slightly dif- 
ferent wavelengths corresponding to these discs introduces 
a modulation with the beat wavelength ~ 2n/Ar. For the 
ring width of 0.5° the wavelength is around 720 which cor- 
responds to the node of the window functions seen in the 
figure. Note also that a larger fraction of the total power for 
rings with a larger radius comes from the lowest multipoles 
in comparison to rings with a smaller radius. Thus, the for- 
mer are more sensitive to the large angular scale power of 
the CMB maps. This dependence will be clearly seen in the 
different tests presented in the paper. 

Finally, it is worth pointing out that one can consider 
a more general estimator than that given by equations 
or |3j). The terms in the sum can be weighted with coeffi- 
cients Wi constrained by additional conditions imposed on 
the estimator. In the case of an optimal estimator, in the 
sense of minimisation of its variance, the weights are given 
by Wi = l/(S + a 2 ) 2 where S = 1/(4tt) + 1) Ci Bj de- 

notes the average CMB variance for the sky observed with an 
azimuthally symmetric beam characterised by the window 
function Bg. Nevertheless, for a high signal-to- noise ratio, 
the variance of the optimal estimator is almost the same 
as for the suboptimal one. This is the case for the 7-year 
WMAP data. While the CMB variance for the V-band map 
is about 7800 /xK 2 , the median noise variance in a pixel is 
about 3300 /xK 2 . Indeed, a comparison between the optimal 
and suboptimal estimator for the WMAP data and found 
only negligible differences. Furthermore, the optimal estima- 
tor requires additional computing time since two additional 
convolutions are performed. For these reasons, we decided 
to used the suboptimal estimator in our analysis. 



4 STATISTICS 

Here we specify the two statistics used to test the signifi- 
cance of the observed variance profiles - the x 2 statistic and 
the number of rings with extreme variance values. 



4.1 x 2 statistic 

In order to test the significance of the measured variances 
for concentric rings centred on a given pixel z, we used the 
X 2 statistic defined as 



Xi 



^ (Var^ r - fjL ijr ) (Var i>r / - jiiy) 



(5) 



where jjLi, r and C^ rr / are the mean variance and covariance 
matrix of the variances for a given pixel z, respectively. Note 
that the application of a mask results in a dependence of 
the mean and covariance matrix values on the position on 
the sky of the ring. Consequently, the computation of the 
Xi statistic becomes complicated. To store all of the covari- 
ance matrices C^ rr / requires of order N p i x N 2 ins bytes of disk 
space. For maps with Af p i x pixels on the sky corresponding 
to a resolution parameter N s id e = 512, and to a number of 
bins A^bins = 180, corresponding to rings of width 0.5° with 
a radius in the range to 90 degrees, this corresponds to 
approximately 400 GB of disk space. Thus to make these 



© 0000 RAS, MNRAS 000, 000-000 



4 P. Bielewicz, B. D. Wandelt, A. J. Banday 




computations feasible, after convolution of the original map 
with the ring window function in spherical harmonic space 
we construct a variance map at N s id e = 128. This is sufficient 
resolution to probe the variance distribution in the sky for 
rings of width 0.5°. Furthermore, we divided the sum over 
bins of ring radii into three intervals: [0°, 30° [, [30°, 60° [ and 
[60°, 90° [. The required disk space is then reduced to around 
2.8 GB. 

The mean and covariance matrices were estimated on 
the basis of 1000 Monte Carlo (MC) simulations of the 7- year 
WMAP maps for the standard A CDM cosmological model 
(see Larson et al.|20l"Tj Table 3) . Another set of simulations 
allowed the estimation of the significance of the observed 
values of the Xi statistic. The Xi statistic was not computed 
for any ring for which more than 50% of the component 
pixels were masked by the WMAP KQ75y7 mask. Such rings 
are delineated by masking the corresponding pixel defining 
the centre of the ring. Masks for the xf maps determined in 
this way are shown in the figures discussed below. 



4.2 Number of extreme variance rings 

Another statistic which we used was the number of rings cen- 
tered on a given pixel with extremely low or high variance. 
Extreme values were defined to be those deviating from the 
mean by more than one and a half standard deviations. To 
estimate the significance of the observed number of extreme 
rings we use MC simulations of the CMB maps and p- values 
for the null hypothesis i.e. the probability of getting a larger 
number of extreme rings for the standard ACDM cosmolog- 
ical model. As was the case for the x 2 statistic, we divided 
the ring radii into three intervals: [0°,30°[, [30°,60°[ and 
[60°,90°[. If concentric rings of different radii centred at a 
given point are independent, then the p- value for the num- 
ber of extreme variance rings n can be estimated from the 
cumulative binomial distribution 



where p is the probability of getting an extreme variance for 
a given ring and N\>i ns = 60 for one of the three radius inter- 
vals of rings with a width of 0.5°. We checked that this ap- 
proximation of the p- value is sufficiently valid for the radius 
interval [0°, 30° [. However, it is not valid for the other inter- 
vals because of strong correlations between rings of larger 
radii. 



5 RESULTS 

Before presenting results of the analysis for the full sky, 
we would like to first address the properties of the sets 
of concentric rings centred on two points which have re- 
cently drawn special attention from the CMB community. 
This interest was based on claims related to the statisti- 
cal significance of the variances measured on these rings in 
the context of the conformal cyclic cosmology scenario pro- 
posed by [Penrose] (2009). The variance profiles for the rings 
centred at the Galactic coordinates (l,b) = (105°, 37°) and 
(l,b) — (252°, —31°) for ring radii from zero to 90 degrees 
are shown in Fig. [2] Comparing the observed profiles with 
simulations based on the standard ACDM model we do not 
identify any significant number of unusually extreme values 
of the variance in good agreement with the conclusions in 



P(n) 



E 

k>n 



A^bins 

k 



v (i -p) 



N hin 



(6) 



Hajian|(|2011| ); |Moss, Scott, fe Zibin] fl2011| ); [Wehus fe Erik- 
|sen| ( |2011 ). However, since our analysis is not limited only 
to a range of radii from to 20 degrees as in the previous 
works, we can extend this conclusion to the full range of ring 
radii from to 90 degrees. 

Nevertheless, for the point (l,b) — (105°, 37°) there is 
evidence of systematically smaller variances for concentric 
rings with radii greater than 10°. Since such rings are more 
sensitive to larger angular scales, we consider that the lower 
variance might be related to the large angular scale anoma- 
lies observed in the WMAP data, and in particular to the 
hemispherical asymmetry in the power spectrum (Eriksen 
et al.|2004| ). Specifically, the point (/, b) = (105°, 37°) lies in 
the northern ecliptic hemisphere which exhibits less power 
than the southern hemisphere. To test this hypothesis, we 
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Figure 2. Variance for rings centred at the Galactic coordinates (I, b) = (105°, 37°) (left and middle) and (7,6) = (252°, —31°) (right) 
with radii in the range [0°,90°[. The red and green lines denote the variance for the V and W-band maps, respectively, masked with 
the KQ75y7 cut. The black lines denote the mean variance estimated from 1000 simulations. The dashed blue lines correspond to the 
mean shifted up and down by one and a half standard deviations of the variance. The shaded dark and light grey regions indicate the 
68% and 95% confidence regions, respectively. In the middle figure, the mean and confidence regions were estimated using simulations 
with multipoles in the range [0, 7] replaced by the corresponding best-fitted multipoles from the foreground corrected V-band WMAP 
map. The dash-dotted line in this figure denotes the mean for simulations without replacing the lowest order multipoles i.e. the same 
mean as in the left figure. Note that the confidence regions are narrower than for the simulations which do not replace the lowest order 
multipoles. This should not be surprising because the cosmic variance for simulations with fixed low order multipoles must be smaller. 



performed additional simulations in which the lowest mul- 
tipoles (£ < 7) were replaced by the corresponding best- 
fit values determined from the foreground corrected V-band 
data for the KQ75y7 sky coverage. Fig. [2] shows that such 
simulations are more consistent with the data. Replacing 
fewer multipoles did not reconcile the anomalous behaviour 
fully, and certainly substituting only the multipoles up to 
order £ = 3 is not sufficient, suggesting that the large scale 
anomalies are not limited only to the apparently aligned 
quadrupole and octopole. 

It is important to note that that a systematically 
smaller variance for rings with larger radii is not consis- 
tent with the conformal cyclic cosmology scenario. In this 
model, low variance rings are generated by the collisions of 
massive black holes in the previous cycle of the Universe. For 
each of the collisions, one should observe a low variance ring 
centred at the position of a given black hole. As the black 
hole undergoes a series of collisions during its lifetime, we 
should observe a series of concentric low variance rings with 
different radii. Earlier collisions will generate rings with a 
larger radius while later ones result in rings with a smaller 
radius. However, for a significant fraction of its lifetime, a 
black hole will not modify the statistical properties of the 
CMB sky and rings not related to the collisions will show 
typical values of the variance. The fact that we find only 
a few rings with variance bigger than the mean for radii 
larger than 10° is not consistent with these predictions. As 
we showed above, the systematically lower variance is much 
better explained by the anomalous distribution of the low 
order multipoles. 

In Fig. [3] we show examples of the variance maps for a 
ring radius of 5 degrees. In appearance, the maps resemble 
a collection of overlapping rings of different amplitude. The 
high amplitude rings are due to pixels with lower or higher 
values in the original map which also form part of the rings 
with centres separated by 5°. It is clear that the variance 
maps are very similar for both V- and W-band maps. Since 
we also did not find significant differences between the re- 
sults for the two bands in the statistics computed below, 



in the reminder of the paper we show results only for the 
V-band map. 

5.1 x 2 statistic 

The x 2_ma P s f° r three radii intervals are shown in Fig. Ul 
The grey regions correspond to pixels for which the % 
statistic was not computed since at least one of the con- 
centric rings centered on a given pixel was masked by more 
than 50%. It is worth noting that in the case of the inter- 
vals [30°, 60° [ and [60°, 90° [ the x 2 statistic is computed for 
those rings centred on some pixels that are located inside 
the KQ75y7 mask. 

There are several regions on the sky with high values 
of x 2 f° r the ring radii in interval [0°,30°[. However, it is 
rather unusual that there is a lack of such regions for the 
X 2 maps corresponding to rings with radii in the intervals 
[30°,60°[ and [60°,90°[. 

To test the significance of the features observed in the 
X 2 maps, we used 1000 MC simulations of the WMAP maps. 
The distributions of the x 2 maps are shown in Fig. [5] If the 
variances for the different rings were uncorrelated, then the 
number of degrees of freedom would be equal to the num- 
ber of bins used in computation of the x 2 i.e. 60 for each 
of the analysed ring radius intervals. However, due to cor- 
relations, the effective number of degrees of freedom will be 
smaller than 60, and this is apparent from the maximum of 
the observed distribution. The distribution for the data lies 
on the border of the 95% confidence region for rings with 
radii in the intervals [30°,60°[ and [60°,90°[. For the inter- 
val [0°,30°[, the data is more consistent with the assumed 
ACDM model. Since rings with larger radii are more sensi- 
tive to large angular scales, it is plausible that the observed 
deviation between data and theory is caused by large an- 
gular scale anomalies. To test this hypothesis, we repeated 
the V-band analysis for the data from which the the best-fit 
multipoles in the range £ £ [0, 7] were removed. Fig.[6]shoes 
that this significantly improves the consistency of the data 
with the model. 
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Figure 3. An example of a variance map for rings with a radius of 5°. The left figure shows the variance for the cross-correlation of 
the VI and V2 DAs computed after application of the KQ75y7 mask. The right figure shows the analogous map determined from the 
cross-correlation of two pairs of W-band DAs - (Wl+W2)/2 and (W3+W4)/2. The grey region corresponds to the pixels for which the 
variance was not computed due to the mask. 



VI x V2 DA, radius range [0,30[ degrees VI x V2 DA, radius range [30,60[ degrees VI x V2 DA, radius range [60,90[ degrees 
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Figure 4. The xf map for the three intervals of the ring radii: [0°,30°[ (left column) , [30°,60°[ (middle column) and [60°,90°[ (right 
column). The maps correspond to the cross-correlation of the VI and V2 DAs masked by the KQ75y7 mask. The black circles are centred 
at two points analysed in Fig. [2] The Grey region corresponds to pixels for which the Xi was n °t computed due to the applied mask. 



Correlations in the variance between rings of differ- 
ent radius are shown in Fig. [7] The correlation matrices 
correspond to rings centred on the Galactic coordinates 
(/,&) = (105°, 37°) for the three ring radius intervals and 
are normalised by the diagonal terms of the covariance ma- 
trix used in the computation of the xl (Eqn.jHJ). Rings with 
a larger radius are strongly correlated. This is seen particu- 
larly well for the intervals [30°, 60° [ and [60°, 90° [. However, 
removing from the map the best-fit multipoles in the range 
i £ [0,7] significantly decorrelates the rings. 

5.2 Number of extreme variance rings 

Fig. [8] shows the probabilities for getting a larger number 
of extreme rings centered at a given point for the V-band 
maps. Regions with unusually numerous extreme variance 
rings, such that none of the 5000 MC simulations show a 
larger number of extreme rings, are denoted by the navy 
blue colour. The colour range from light blue to red denotes 
regions with a typical or smaller than typical number of 
extreme variance rings for the ACDM model. 

In the case of the interval [0°,30°[, we can notice that 
over a large part of the sky, and especially in the northern 
Galactic hemisphere, the number of rings with extremely 
high variance is slightly smaller than for most of simulations. 
However, one can find also regions, mostly in the south- 
ern Ecliptic hemisphere, where the number of high variance 
rings is much larger than for the simulations. As one would 



expect, these regions are clearly anti-correlated with the 
parts of the sky corresponding to rings with extremely low 
variance seen in the lower row of the figure i.e. regions with a 
higher number of extremely high variance rings correspond 
to regions with lower number of extremely low variance rings 
and vice- versa. 

One should notice that the point (/,&) = (105°, 37°), 
corresponding to the center of a set of concentric rings with 
systematically smaller variances, lies in one of the region 
close to the north Ecliptic pole with many extremely low 
variance rings. Any other point in this region shows a similar 
number of unusually low variance rings. Such a clustering of 
rings with low variance is not consistent with the conformal 
cyclic cosmology scenario. Since the positions correspond- 
ing to the massive black holes which generate low variance 
rings through collisions are distributed randomly on the sky, 
the centers of the rings should also be distributed randomly. 
Clustering of the rings in relatively large regions rather in- 
dicates effects related to large angular scale features in the 
CMB maps. 

The anti-correlation between high variance and low 
variance ring regions is seen also for the intervals [30°,60°[ 
and [60°,90°[. However, in these cases the p- values for the 
number of extremely high variance rings are smaller than for 
the interval with radii [0°, 30° [.This is a consequence of the 
stronger correlations between rings of different radius, as the 
variance for these intervals is more sensitive to the large an- 
gular scales. One can also notice fewer regions with small p- 
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VI x V2 DA, radius range [0,30[ deg VI x V2 DA, radius range [30,60[ deg VI x V2 DA radius range [60,90[ deg 




Figure 5. The distribution of the xf map for the three intervals of the ring radii: [0°,30°[ (left) , [30°,60°[ (middle) and [60°,90°[ 
(right). The red lines denote the distribution for the WMAP maps masked with the KQ75y7 mask and the black lines indicate the mean 
of the distribution derived from 1000 MC simulations of signal plus noise. The shaded dark and light grey regions indicate the 68% and 
95% confidence regions, respectively. The distributions correspond to the cross-correlation of the VI and V2 DA maps. 



VI x V2 DA radius range [0,30[ deg VI x V2 DA radius range [30,60[ deg VI x V2 DA radius range [60,90[ deg 




40 60 80 100 120 140 40 60 80 100 120 140 40 60 80 100 120 140 



Figure 6. The distribution of the xf map after removal of the best-fit multipoles in the range t £ [0, 7] from the VI and V2 DAs. Plotted 
lines as in Fig. [5] 




Figure 7. Correlation matrices for the variance on concentric rings centred at the Galactic coordinates (l,b) = (105°, 37°) and ring 
radius intervals [0°,30°[ (left column) , [30°,60°[ (middle column) and [60°,90°[ (right column). The bottom row corresponds to the 
matrices derived after removing the best-fit multipoles in the range i G [0, 7] from the maps. The matrices were estimated using 1000 
MC simulations. 
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values than for the interval [0°, 30° [. Conversely, over a large 
part of the sky and especially for the interval [30°, 60° [ and 
the northern Ecliptic hemisphere, the number of extremely 
low variance rings is significantly larger than for simulations. 
In the case of the interval [60°,90°[ the low variance rings 
are centered mostly in the vicinity of the Galactic center. 
This disparity in size of regions with extremely high and 
low variance rings is related to the deficit of power at large 
angular scales in comparison to the best-fit ACDM model. 

For comparison we show in Fig. [9] the p- value maps for 
simulated VI and V2 DA maps with the amplitude of the 
quadrupole close to the average predicted for the concor- 
dance ACDM model. As we see, contrary to the WMAP 
data, in this case regions with small p- value are much larger 
for extremely high variance rings and much smaller for ex- 
tremely low variance rings. This is seen very well for the 
radius intervals [30°,60°[ and [60°,90°[. This comparison 
confirms again the deficit of power of the WMAP maps at 
large angular scales. It also shows that relatively large re- 
gions with very small or large probability of getting larger 
number of extreme rings are generic. Thus, the p- value maps 
are rather not useful for estimation of how significant is this 
deficit. However, they can be used to find parts of the sky 
with a large number of extreme variance rings that may be 
interesting for the search for signatures of alternative models 
of structure formation such as conformal cyclic cosmology. 

As for the x 2 map, we repeated estimation of the p- 
value maps after removing the best-fit multipoles in the 
range £ £ [0,7] from both simulations and the V-band DAs. 
The results are shown in Fig. |10| Since the variance for rings 
with radii in the interval [0°, 30° [ is the least sensitive to the 
large angular scales, removal of the lowest order multipoles 
does not significantly change the p-value maps for this in- 
terval. The most noticeable change is the reduction of the 
region of the low variance rings close to the north Ecliptic 
pole. More significant changes are noticeable for the other 
two intervals. Since the correlations between rings with dif- 
ferent radii after removing of the lowest order multipoles are 
smaller, the distribution of the p- values looks more similar 
to the distribution for the interval [0°,30°[. In the case of 
the high variance rings, we can also notice extended regions 
with lower p- values in the southern Galactic hemisphere. On 
the other hand, regions corresponding to the extremely low 
variance rings are smaller. 

Finally, we undertook another test of dependence of the 
results on the low order multipoles. We replaced the best-fit 
multipoles in the range £ £ [0,7] by the same multipoles ro- 
tated around the z-axis by 120°. The results for such maps 
are shown in Fig. After rotation, regions with extremely 
low variance rings, especially for the intervals [30°, 60° [ and 
[60°,90°[, are significantly smaller. However, for the high 
variance rings, the changes are less significant. This test 
would appear to again confirm that, to a large extent, the 
observed extreme variances of the rings are caused by the 
lowest order multipoles. 



6 CONCLUSIONS 

We have presented a fast method for the computation of 
variance on azimuthally symmetric patches of the sky using 
a fast convolution approach. This has enabled the search 



for a set of concentric rings of width of 0.5° centred at the 
nodes of a higher resolution grid than would be possible for 
the pixel-based approach used in previous studies. It has 
also facilitated the analysis of rings over the full range of 
radii from to 90 degrees. The technique was applied to the 
7- year WMAP data to search for rings with unusually low 
or high variance. 

Analysis of the variance of rings centred at two points, 
(Z,6) = (105°, 37°) and (/, b) = (252°, -31°), which have 
recently drawn special attention in the context of the con- 
formal cyclic cosmology scenario, showed that for the former 
with ring radius larger than about 10° the variance is sys- 
tematically smaller than the average from MC simulations. 
However, after replacing the £ £ [0,7] multipoles in the simu- 
lations with the best-fit values from the foreground corrected 
V-band data masked with the KQ75y7 cut, the agreement 
between simulations and the observed variances is very good. 
Therefore, we can conclude that observed deficit of power 
in rings with large radii is related to the anomalous distri- 
bution of the lowest order multipoles of the WMAP map. 
Moreover, the p-value maps for the number of extremely 
low variance rings show that one can easily find points with 
as many anomalously low variance rings as for the point 
(l,b) = (105°, 37°) in large part of the northern ecliptic 
hemisphere. Thus, there is nothing unusual in the properties 
of the rings centred at the point (Z, b) = (105°, 37°) given 
the observed large scale anomalies. These results are not 
consistent with predictions of the conformal cyclic cosmol- 
ogy scenario and they can not be used to support it. Whilst 
this model predicts randomly distributed sets of concentric 
rings within which few have lower variance, we observe sets 
of rings with systematically lower variance for larger ring 
radii that are clustered in a few relatively large regions of 
the sky 

A strong dependence of the variance in rings, especially 
with large radius, on the lowest order multipoles is also seen 
for other statistics used in our studies. The \ 2 map cor- 
responding to the radius intervals [30°,60°[ and [60°,90°[ 
takes smaller values than simulations for the best-fit ACDM 
model at a significance level around 95%. Furthermore, in a 
large part of the northern Ecliptic and Galactic hemispheres, 
the probability of getting a larger number of extremely low 
variance rings is very low, while extremely high variance 
regions are much smaller and localised mostly in the south- 
ern Ecliptic hemisphere. All of these anomalies disappear 
after removing from the maps the best-fit multipoles in the 
range £ £ [0, 7]. The disparity in distribution of the low and 
high variance rings is also weaker after rotation of the low- 
est order multipoles around the Galactic poles. Thus, the 
anomalous variances in rings can be traced to the presence 
of anomalous low order multipoles. 

However, as we tested, removal of the quadrupole and 
octopole is insufficient to eliminate completely the anoma- 
lies. Therefore, we conclude that they are related also with 
the higher order multipoles, up to £ = 7. This conclusion is 
consistent with the results of previous studies on the large 
angular scale anomalies observed in the WMAP data (Hou, 



Banday, & Gorski 2009 Hou et al. 2010). Studies of the 



physical origin of these anomalies is beyond the scope of this 
paper. However, they do not seem to be related to residu- 
als of the Galactic thermal dust emission since the results 
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Figure 8. The p-values map, estimated from 5000 MC simulations, for rings with extremely high (upper row) and low (lower row) 
variances for the three intervals of the ring radii: [0°,30°[ (left column) , [30°,60°[ (middle column) and [60°,90°[ (right column). The 
variance map was estimated using the VI and V2 differencing assemblies WMAP maps. The black circles are centred at two points 
analysed in Fig. [2] The grey region corresponds to pixels for which the p-values were not computed due to the applied mask. 
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Figure 9. The same as in Fig. |8| but for simulated VI and V2 DAs. 



for the W-band map are consistent with the results for the 
V-band map. 
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Figure 10. The same as in Fig. ||] but for maps with the best-fit multipoles in the range t G [0, 7] removed. 
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Figure 11. The same as in Fig.[8]but after rotating the best-fitted multipoles in the range t G [0, 7] by 120° around the z-axis. 
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